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Accurate and early detection of plant diseases will facilitate mitigate the 
worldwide losses experienced by the agriculture area. MATLAB image 
processing provides quick and non-destructive means of rust disease 
detection. In this paper, microscopic image data of rust disease of Lentil was 
combined with image processing with depth information and developed a 


machine learning system to detect rust disease at early stage infected with 


fungus Uromyces fabae (Pers) de Bary. A novel feature set was extracted 
Keywords: from the image data using local binary pattern (LBP) and HBBP (Brightness 
Bi-Histogram Equalization) for image enhancement. It was observed that by 
combining these, the accuracy of detection of the diseased plants at 
; microscopic level was significantly improved. In addition, we showed that 
Lentil our novel feature set was capable of identifying rust disease at haustorium 
Local binary pattern stage without spreading of disease. 
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1. INTRODUCTION 

India is an agriculture based country and agricultural productivity directly affects the economical 
growth of a developing country. Plant diseases are the major concern in crop productivity. Automatic 
detection of disease in plants is an emerging area for contemporary day biocomputing science. Image 
processing and pattern recognition techniques have paved the way for accurate and early detection of fungal 
diseases in agriculture [1-3]. MATLAB image processing has been used by various researchers to study the 
effect of disease on a plant [4-7]. 

Among various imaging techniques, digital imaging has been shown to be a robust technique for 
identification of fungal diseased regions on leaves [8, 9]. One amongst the main issue related to digital 
imaging in plants is the late detection of rust disease because at visual stage rust disease could easily spread 
with insects and wind to other fields and couldn’t be controlled at an early stage. So we need better technique 
to close the gap between visualization of disease and objective quantitative results to control rust disease at 
an early stage before spreading to other fields. Microscopic images could be better alternatives to detect rust 
disease at early i.e. haustorium stage [10, 11]. Seiffert et al. (2005) used pattern recognition software for 
automatic detection of hyphal growth of powdery mildew disease of barley i.e. Blumeria graminis f. sp. 
Hordei [12]. 

In this paper, we tend to combine microscopic images with image processing techniques to 
overcome these issues and present a method for automatic detection of rust disease at an early stage using 
machine learning techniques. This work is a step towards creating a novel approach for automatic detection 
of rust disease of pea plants caused by fungus Uromyces fabae (Pers) de Bary before the appearance of 
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visible symptoms or before spreading of disease. For depth estimation, we proposed two complimentary 
approaches Local binary patterns (LBP) to extract a novel feature set and BBHE (Brightness Bi-Histogram 
Equalization) for image enhancement to show that it is capable of automatic detection of leaves poised to be 
suffering from rust throughout the experiment. 


2. RESEARCH METHOD 
2.1. Image acquisition 

An experimental setup was designed and developed to simultaneously acquire microscopic images 
of leaves to distinguish leaves with or without rust disease. The imaging setup consisted of a light microscope 
with the camera (LIECA). The experiment was carried out on 300 images of leaves were collected for 15 
consecutive days and stained with congo red dye to improve the appearance of fungal haustorium. The 
disease symptoms that developed consisted of haustorium that expanded over time in the form of spores and 
eventually caused chlorosis and ready for spreading to other plants as well as fields. Figure la showed the 
haustorium of rust after staining. 


2.2. Selection of an image 
Initially, an image had selected by the user on which the operation was performed. In MATLAB, 
“uigetfile” function was used to browse an image from everywhere in the computer. The syntax used to get 
the image was as follows: 
[File Path]=uigetfile (‘.jpg', ‘Select An Image’); 
Where the get image was stored in a matrix in the form of the file name and corresponding path. 
After selecting an image, a few steps were performed to make it feasible for further procedures. 
a. Read the selected image using “imread” function available in the MATLAB. 
Resize the read image in order to reduce the processing time. 
c. Lastly, show the selected image considered as an original image using “imshow’” function of the 
MATLAB 


2.3. Application of median filter 

Once the image had selected and output had shown, a median filter was applied over the selected 
image. The application of this filter was smoothening of the actual image while preserving the edges. The 
median filter uses a window size where one entry immediately preceding and following each entry. In the 
proposed work, the median filter was applied to the image to smoothening of the image as shown in 
Figure 1b which had obtained after applying median filter onto the original image. 
Allocate the actual size of image i.e. width, height, and layer. 
Assign the size of the window where its width and height is initialized. 
Set for loop for the number of layers from | to total size of the image. 
Initialize the loop for rows from 1 to total size of the image rounded down by 2. 
Initialize the loop for Columns from 1 to total size of the image rounded down by 2. 
e.l. Now initialize the for loop x from 1 to 3 
e.2. Initialize the for loop y from 1 to 3 
Shift the window over the output image one by one and increment the location of window by 1. 
f. End of the for loop y 
g. End of the for loop x 
h 
i. 


ono Ff 


. Evaluate the median filter of the obtained window. 
Demonstrate the acquired smoothened image. 


2.4. Apply BBHE enhancement 

Once the smoothened image had acquired from the previous step, next step was to enhance the 
image. For the enhancement purpose, BBHE technique was used. The BBHE i.e. Brightness Bi-Histogram 
Equalization algorithm preserves the mean brightness of a given image while enhancing the contrast of that 
image. After apply BBHE for the enhancement purpose over the image, the acquired enhanced output image 
has shown in Figure 1c and Figure 1d. 
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Figure 1. (a) Original Image; (b) Smoothened Image using Median Filter; (c) enhanced image after applying 
BBHE; (d) Feature extracted from an image using Local Binary pattern 


Divide the histogram image into two sub parts 
Perform normalization using PDF and CDF 


Combine two histograms and acquired final equalized 
image 


Figure 2. Flowchart of the BBHE algorithm 
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2.5. Algorithm for BBHE 
a. Evaluate the mean of the smoothened image and considered as a mean threshold value. 
b. Generate two histograms i.e. lower histogram and upper histogram for the evaluation of a number of 
count of each color. 
Set a for loop i for number of rows. 
Set a for loop j for number of columns. 
e. Now compare the pixel value of an image with the mean threshold value and decomposed the original 
image into sub-images as lower_histogram and Upper_histogram. 
1. Ifthe mean threshold value is greater than the pixel’s value than it will be lower histogram. 
2. Ifthe mean threshold value is lesser than the pixel’s value than it will be upper histogram. 
f. End of the for loop j 
End of the for loop i 
Calculates the PDF i.e. Probability Density Function of the lower and upper histogram for 
normalization. 
1. Normalize_lower_histo=histo_lower/sum(histo_lower); 
2. Normalize _upper_histo=histo_upper/sum(histo_upper); 
i. Evaluate the Cumulative Density Function for lower and upper range. 
j. Initiate the for loop K from 2 to 255 and estimates the lower and upper histogram’s CDF. 
1. hist_lower_cdf(k)=hist_lower_cdf(k-1) +Normalize_lower_histo (k); 
2. hist_upper_cdf(k)=hist_upper_cdf(k-1) +Normalize_upper_histo (k); 
k. Initiate the loop i from 1 to total number of rows. 
1. Initiate the loop j from 1 to total number of columns. 
m. Define transform function based on the CDF such as: 
1. final_image(i,j)=range_lower(1)+round(((range_lower(2)- 
range_lower(1))*hist_lower_cdf(original_pixel_value+1))); 
2. final_image(i,j)=range_upper(1)+round(((range_upper(2)- 
range_upper(1))*hist_upper_cdf(original_pixel_value+1))); 
n. The final attained histogram equalized image will be expressed as: 
EnhancedImage(:,:,ij)=Equilzed_ BBHE (Ouptut Image (:,:,1j)); 


eo 


= ga 


2.6. Feature extraction using LBP 
After enhancing the quality of the input image, features were extracted from that image. For the 
feature extraction purpose, in the proposed work, Local Binary Patterns operator was used. This operator was 
basically used to describe the texture and shape of a digital image. The process was done by dividing the 
whole image into small regions. Further, these regions were used to extract the features as shown in Figure 
1(d) and flow chart of LBP in Figure 3. The steps followed by the LBP were: 
a. Initialize the size of the LBP label “L”. 
b. Compute the center of the label by dividing it with 2. 
C=round (L/2) 
c. Initialize the size of the window by computing row_max and col_ max. 
row_max=size(Input_Image,1)-L+1; 
col_max=size(Input_Image,2)-L+1; 
d. Compare each pixel of the cell to its neighbors such as Left-top, Left-middle, left bottom, right-top, 
right-middle, right-bottom etc. Follow the pixels along a circle such as clockwise or counterclockwise. 
e. Compare the center pixel’s value, if greater than “0” write “1”. This process will end up with 8-digit 
binary number. 
f. | Convert the generated binary number into decimal using below equation. 
LBP_Image(i,j)=(Center,Left) +(Left,Left)*2 +(Left,Center)*4 +(Left,1)*8 +(Center,1)*16 +(1,1)*32 
+(1,Center)*64 +(1,Left)*128; 
g. The attained pixel is the LBP output image. 
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Initialize the size of the window 
Evaluate the center of the window 
Compare each pixel of the cell with its neighbor 


If neighbor’s pixel 
value >0 


Write 1 on the position 
Covert the attained binary number into decimal 


Figure 3. Flowchart of the LBP feature extraction technique 


2.7. Evaluate the color value to estimates the intensity value of the input image 
In this step, the intensity of the image was evaluated. Thus, the steps followed for the same were: 
a. Evaluate the size of the image. 
b. Ifthe numbers of layers are three then extract each layer from an image i.e. red, green and blue. 
r=image (:,:,1); % red color 
g=image (:,:,2); % green color 
b=image (:,:,3); % blue color 
c. Calculate the mean of each layer individually: 
red=mean(mean(r)); 
green=mean(mean(g)); 
blue=mean(mean(b)); 
d. Nowevaluate the mean of all the layers: 
main_value=(red+green+blue)/3; 
e. Resultant value will be the color value. 


2.8. Identify the shape value of the image 
This step provided the shape value of an image in the form of numeric value. The higher the Value, 
the chances of having circle shape were also high. Algorithm to evaluate the shape was as follows: 
a. Take the RGB image. 
b. Convert the RGB image into gray format. 
c. Evaluate the threshold value of an image. 
Threshold=graythresh(image); 
Convert the threshold image into black and white image. 
e. Remove the noise by eliminating the entire object that contains fewer than 30 pixels. 
f. Find the corresponding boundaries 
f.1. Display the label matrix and draw each boundary using the for loop k. 
boundary=B {k}; 
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g. Determine the round objects by evaluating the area and center using regionprops function. 

1. Initiate a for loop from 1 to total length of the boundaries. 
Obtain boundary coordinates. 

3. Compute the perimeter of an object: 
delta_square=differences(boundary).”2; 
perimeter=sum(square root(sum(delta_sq,2))); 

4. Obtain the total area calculated. 

5. Compute the roundness of the metric. 
Metric=4*pi*area/perimeter’2; 

6. End of the for loop 

h. Acquired the shape value of an image. 


3. RESULTS AND ANALYSIS 

In this work, image processing of rust disease was performed as shown in Figure 4. At the time of 
acquisition, the image may be corrupted due to noise that degrades the actual quality of the image. So, in 
order to extract the features, it should be enhanced firstly. In the proposed work, initially, pre-processing 
steps have followed where the image is filtered and enhanced and then features are extracted from the 
enhanced image as shown in flowchart of the proposed methodology. The output image that has acquired 
after apply BBHE algorithm over the smoothened image. From the result, it has been cleared that the output 
image is more enhanced and bright in comparison with the image which has attained in the previous step as 
shown in Figure 4. After feature extraction color and shape value of the input image were found 191.5965 
and 0.7213, respectively. Previously, automatic detection methods for rust disease of the Wheat crop were 
reported [13, 14]. Image enhancement and feature extraction were reported by many scientists for the fungal 
disease of crops [1, 2, 15, 16]. 
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Figure 4. Framework of the proposed model 


4. CONCLUSIONS 

This paper exhibited the utilization of microscopic images to detect rust disease on lentil leaves. The 
different lentil varieties did not influence the algorithm’s ability to recognize rust disease on leaves. BBHE 
techniques produced the best image contrast enhancement, without unfortunate antiquities and maintain input 
mean brightness of images. This investigation endeavored to check the achievability of an image processing 
technique for rust disease detection and the outcomes introduced are promising. Future work is aimed to 
carry out a similar experiment to detect rust of other crops. One of the fundamental difficulties looked in this 
investigation was staining protocol which causes a commotion in images. The primary target will be to 
modify and improve the algorithm to develop a robust algorithm that can make exact recognition of rust 
disease in field conditions. 
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